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The diffuse Galactic 7-ray emission is produced by cosmic rays (CRs) interacting with the inter- 
stellar gas and radiation field. Measurements by the Energetic Gamma-Ray Experiment Telescope 
(EGRET) instrument on the Compton Gamma-Ray Observatory indicated excess 7-ray emission 
> 1 GeV relative to diffuse Galactic 7-ray emission models consistent with directly measured CR 
spectra (the so-called "EGRET GeV excess" ) . The excess emission was observed in all directions on 
the sky, and a variety of explanations have been proposed, including beyond-the-Standard-Model 
scenarios like annihilating or decaying dark matter. The Large Area Telescope (LAT) instrument 
on the Fermi Gamma-ray Space Telescope has measured the diffuse 7-ray emission with improved 
sensitivity and resolution compared to EGRET. We report on LAT measurements of the diffuse 
7-ray emission for energies 100 MeV to 10 GeV and Galactic latitudes 10° < |6| < 20°. The LAT 
spectrum for this region of the sky is well reproduced by a diffuse Galactic 7-ray emission model 
that is consistent with local CR spectra and inconsistent with the EGRET GeV excess. 

PACS numbers: 95.30.Cq,95.55.Ka,95.85.Pw,96.50.sb,98.70.Sa 
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Introduction: The diffuse 7-ray emission, both Galac- 
tic and extragalactic, is of significant interest for as- 
trophysics, particle physics, and cosmology. The dif- 
fuse Galactic emission (DGE) is produced by interac- 
tions of cosmic rays (CRs), mainly protons and elec- 
trons, with the interstellar gas (via 7r°-production and 
bremsstrahlung) and radiation field (via inverse Compton 
[IC] scattering) It is a direct probe of CR fluxes in 

distant locations, and may contain signatures of physics 



beyond the Standard Model, such as dark matter annihi- 
lation or decay. The DGE is a foreground for point-source 
detection and hence influences the determination of the 
source positions and fluxes. It is also a foreground for the 
much fainter extragalactic component, which is the sum 
of contributions from unresolved sources and truly diffuse 
emission, including any signatures of large scale structure 
formation, emission produced by ultra-high-energy CRs 
interacting with relic photons, and many other processes 
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(e.g., Q and references therein). Therefore, understand- 
ing the DGE is a necessary first step in all such studies. 

The excess diffuse emission > 1 GcV in the Ener- 
getic Gamma-Ray Experiment Telescope (EGRET) data 
[J] relative to that expected from DGE models consis- 
tent with the directly measured CR nucleon and elec- 
tron spectra 0, @] led to the proposal that this emis- 
sion was the long-awaited signature of dark matter anni- 
hilation Q. More conventional interpretations included 
variations of CR spectra in the Galaxy [!, [H, 0] , contribu- 
tions by unresolved point sources [8|, and instrumental 
effects [IIIEkI. 

A model of the DGE depends on the CR spectra 
throughout the Galaxy as well as the distribution of the 
target gas and interstellar radiation field (1SRF). Start- 
ing from the distribution of CR sources and particle in- 
jection spectra, the distribution of CRs throughout the 
Galaxy is determined taking into account relevant energy 
losses and gains, then the CR distributions are folded 
with the target distributions to calculate the DGE [e.g., 
111 ]. Defining the inputs and calculating the models are 
not trivial tasks and involve analysis of data from a broad 
range of astronomical and astroparticle instruments [13] • 

The Fermi Large Area Telescope (LAT) was launched 
on June 11, 2008. It is over an order of magnitude more 
sensitive than its predecessor, EGRET, with a more sta- 
ble response due to the lack of consumables. The LAT 
data permit more detailed studies of the DGE than have 
been possible ever before. 

In this paper, analysis and results for the DGE are 
shown for the Galactic mid-latitude range 10° < \b\ < 20° 
measured by the LAT in the first 5 months of the science 
phase of the mission. This region was chosen for initial 
study since it maximises the fraction of signal from DGE 
produced within several kpc of the Sun and hence un- 
certainties associated with CR propagation, knowledge 
of the gas distribution, etc., should be minimised. The 
calculation of the DGE at lower Galactic latitudes re- 
quires CR fluxes throughout the whole Galaxy and thus 
is model dependent, while the emission at higher lati- 
tudes is more affected by contamination from charged 
particles misclassified as photons and uncertainties in the 
model used to estimate the DGE. The diffuse emission at 
lower and higher Galactic latitudes will be addressed in 
subsequent LAT papers. 

LAT Data Selection and Analysis: The LAT is a 
pair-conversion telescope with a precision tracker and 
calorimeter, each consisting of a 4 x 4 array of 16 mod- 
ules, a segmented anti-coincidence detector (ACD) that 
covers the tracker array, and a programmable trigger and 
data acquisition system. Full details of the instrument, 
onboard and ground data processing, and other mission- 
oriented support are given in [131 ]. 

The data selection used in this paper is made using the 
standard LA T g round processing and background rejec- 
tion scheme 13j. This consists of two basic parts: first 



a simple accept-or-reject selection (prefiltering) followed 
by a classification tree (CT) [3] based determination of 
the relative probability of being background or signal. 
The prefiltering phase screens particles entering the LAT 
for their charge neutrality using the tracker and ACD. 
The direction reconstruction software extrapolates par- 
ticle trajectories found in the tracker back to the scin- 
tillation tiles of the ACD, and we accept only events in 
which the intersected tiles show no significant signal. In 
addition, the prefiltering phase includes considerations 
of the shape of the calorimeter shower energy deposition 
and how well the found tracks project into the energy 
centroid. The overall background rejection of the pre- 
filtering phase is 10 3 — 10 4 depending on energy, yielding 
an efficiency > 90% for 7-rays that convert into electron- 
positron pairs in the LAT. 

Classification trees, which afford an efficient and sta- 
tistically robust method for distinguishing signal from 
noise, are used to reduce backgrounds further. Using 
quantitites defined from ACD, tracker, and calorimeter 
data, the CTs are trained on Monte Carlo simulated data 
which have passed the prefilter described above. Multiple 
CTs are built to make the procedure robust against sta- 
tistical fluctuations during the training procedure. The 
result from averaging the output from these CTs is the 
probability for an event to be a photon or background. 
This final selection parameter allows the signal purity to 
be set according to the needs of the analysis. For the 
analysis of diffuse emission, the cut on the CT generated 
probability is set such that the Monte Carlo prediction 
of the orbit-averaged background rate is ~ 0.1 Hz inte- 
grated over the full instrument acceptance > 100 MeV. 
This yields a 7-ray efficiency > 80%, and the residual 
background is at a level where the majority of the con- 
tamination arises from irreducible sources such as 7-rays 
produced by CR interactions in the passive material out- 
side the ACD, e.g., the thermal blanket and microme- 
teroid shield of the LAT (see Fig. 13 in [Tj|). The events 
corresponding to the above criterion are termed "Diffuse" 
class and are the standard low-background event selec- 
tion. 

The analysis presented here uses post-launch instru- 
ment response functions (IRFs). These take into account 
pile-up and accidental coincidence effects in the detec- 
tor subsystems that are not considered in the definition 
of the pre-launch IRFs. Cosmic rays, primarily protons, 
pass through the LAT at a high rate and sufficiently near 
coincidences with 7-rays leave residual signals that can 
result in 7-rays being misclassified, particularly at en- 
ergies < 300 MeV. The post-launch IRFs were derived 
using LAT events read from a special trigger that pro- 
duces periodic detector readouts, irrespective of the sig- 
nals present, as a background overlay on the standard 
simulations of 7-rays and provide an accurate accounting 
for the instrumental pile-up and accidental coincidence 
effects. The on-axis effective area for the event selection 
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used in this paper is ~ 7000 cm 2 at 1 GeV and is energy 
dependent; this is approximately 10% lower at 1 GeV 
than the pre-launch effective area corresponding to the 
same event selection. The systematic uncertainties of the 
effective area, evaluated by comparing the efficiencies of 
analysis cuts for data and simulation of observations of 
Vela, are also energy dependent: 10% below 100 MeV, 
decreasing to 5% at 560 MeV, and increasing to 20% at 
10 GeV and above. The point spread function (PSF) and 
energy resolution are as described in [l3j ]. 

The LAT nominally operates in a scanning mode that 
covers the whole sky every two orbits (i.e., 3 hrs). We 
use data taken in this mode from the commencement 
of scientific operations in mid-August 2008 to the end 
of December 2008. The data were prepared using the 
LAT Science Tools package, which is available from the 
Fermi Science Support Center d5]. Events satisfying 
the Diffuse class selection and coming from zenith an- 
gles < 105° (to greatly reduce the contribution by Earth 
albedo 7-rays) were used. To further reduce the effect of 
Earth albedo backgrounds, the time intervals when the 
Earth was appreciably within the field of view (specifi- 
cally, when the centre of the field of view was more than 
47° from the zenith) were excluded from this analysis. 
This leaves 9.83 Ms of total livetime in the data set. 
The energy-dependent exposure was calculated using the 
IRFs described above. 

The photon counts and exposure were further pro- 
cessed using the GaRDiAn package, part of a suite of 
tools we have developed to analyse the DGE where the 
analysis approach is described in [16J, with more details 
to be given in a subsequent publication. Gamma-ray 
skymaps were generated using a HEALPix [13] scheme 
at order 7 (i.e., ~ 0.2° resolution) with 5 bins per decade 
in energy from 100 MeV to 10 GeV. For each energy bin 
the intensity was obtained by dividing the in-bin counts 
by the spectrally-weighted exposure over the bin. We 
used two methods for the spectral weighting: a power 
law with index —2 and the spectral shape of the assumed 
DGE model (described below). With the energy binning 
used in this paper the differences in the derived intensi- 
ties were < 1% between these two weighting schemes. 

Figure[T]shows the LAT data averaged over all Galactic 
longitudes and the latitude range 10° < \b\ < 20°. The 
hatched band surrounding the LAT data indicates the 
systematic uncertainty in the measurement due to the 
uncertainty in the effective area described above. Also 
shown are the EGRET data for the same region of sky 
derived from count maps and exposures available via the 
CGRO Science Support Center fl 811 and processed follow- 
ing the procedure described in [11[ and we have included 
the standard systematic uncertainty of 13% jl9| . For 
both data sets the contribution by point sources has not 
been subtracted. The LAT-measured spectrum is sig- 
nificantly softer than the EGRET measurement with an 
integrated intensity J L at(> 1 GeV) = 2.35±0.01xl0~ 6 



cm 2 s 1 sr 1 compared to the EGRET integrated in- 
tensity J E gret(> 1 GeV) = 3.16 ±0.05 x 10~ 6 cm" 2 s" 1 
sr -1 where the errors are statistical only. Not included 
in the figure is the systematic uncertainty in the energy 
scale, which is conservatively estimated from comparison 
between Monte Carlo and beam test data as < 5% for 
100 MeV to 1 GeV, and < 7% above 1 GeV where it is 
believed that if any bias is present energies are overesti- 
mated. Taking the uncertainty on the energy scale into 
account, the LAT spectrum could be softer, increasing 
the discrepancy with the EGRET spectrum further. 
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FIG. 1: Diffuse emission intensity averaged over all Galactic 
longitudes for latitude range 10° < \b\ < 20°. Data points: 
LAT, red dots; EGRET, blue crosses. Systematic uncertain- 
ties: LAT, red; EGRET, blue. 

Figure [3 compares the LAT spectrum shown in Fig.[T] 
with the spectra of an a priori DGE model, and a point- 
source contribution and unidentified background (UIB) 
component derived from fitting the LAT data that are 
described below. The DGE model is an updated version 
of the "conventional" model from GALPROP [nj]. Ma- 
jor improvements include use of the formalism and corre- 
sponding code for pion production in pp-interactions by 
[201 ] , a complete recalculation of the ISRF 21 1 , updated 
gas maps, and an improved line-of-sight integration rou- 
tine. However, it is still an a priori model that is based 
on local cosmic-ray data, and does not use 7-ray data. 
Table U summarises the numerical values by energy bin 
for the different components shown in Fig. [2j 

The source and UIB components were obtained by fit- 
ting the LAT data using GaRDiAn with the DGE model 
held constant. Point source locations were taken from 
the 3 month Fermi LAT source list down to sources with 
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FIG. 2: LAT data with model, source, and UIB compo- 
nents for sky region in Fig. [Tj Model (lines): 7r°-decay, 
red; bremsstrahlung, magenta; IC, green. Shaded/hatched 
regions: UIB, grey/solid; source, blue/hatched; total (model 
+ UIB + source), black/hatched. 



5-<7 significance. Due to the limited statistics of all but 
the very brightest sources, we used 3 bins per energy 
decade in the fitting procedure. Source positions were 
fixed but the spectra were fit using one free parameter 
for the source flux per energy bin. The UIB component 
was determined by fitting the data and sources over all 
Galactic longitudes for the high-latitude region |6| > 30° 
for the full LAT energy range shown in the figure. Using 
this high-latitude region minimises the effect of contam- 
ination by the bright Galactic ridge which can be signif- 
icant even up to ~ 10° from the plane due to the long 
tails of the PSF at low energies. 

To determine the uncertainty of the source and UIB 
components, we modified the effective area to the ex- 
tremes of its systematic uncertainty defined before and 
refitted the data. Since the DGE model components do 
not vary in the fit, the absolute change in intensity caused 
by the modification to the effective area propagates di- 
rectly to the source and UIB components. The system- 
atic uncertainty on these components is energy depen- 
dent and due to several effects. 

For energies > 10 GeV the PSF is - 0.2° (68% contain- 
ment) and the sources are well-localised spatially. Since 
the model is fixed and the sky maps are sparser at high 
latitudes for the data taking period in this paper, the 
UIB component absorbs almost all of the intensity from 
the modification to the effective area. At low energies 
the PSF is wider, 3.5° (68% containment) at 100 MeV 



TABLE I: 

Energy a 



LAT data and components: 



LAT 



be 



Model 



b c d 



10° < |6| < 20°. 

UIB 6e/ Source 66 



100-158 59.8 ±0.3 26.0 11.0 6.4 8.6 21.0 ±0.1 7.4 ± 0.4 

158-251 65.0 ±0.3 33.5 18.2 7.3 8.0 20.5 ± 0.1 5.4 ± 0.1 

251-398 67.1 ±0.3 38.2 23.2 7.6 7.4 18.7 ± 0.1 4.7 ± 0.1 

398-631 64.5 ± 0.3 38.9 25.3 7.0 6.6 15.4 ± 0.1 3.8 ± 0.1 

631-1000 60.8 ±0.3 37.3 25.7 5.7 5.9 12.9 ± 0.1 3.2 ± 0.1 

1000-1585 55.1 ±0.4 32.8 23.3 4.4 5.1 11.6 ± 0.1 3.0 ± 0.1 

1585-2512 46.3 ±0.4 26.5 19.0 3.1 4.4 10.7 ±0.1 2.7 ±0.1 

2512-3981 37.0 ± 0.5 20.2 14.4 2.0 3.8 9.2 ± 0.1 2.3 ± 0.1 

3981-6310 29.9 ±0.5 14.9 10.5 1.2 3.2 8.5 ±0.1 1.9 ±0.1 

6310-10000 20.7 ±0.5 10.9 7.5 0.7 2.7 6.8 ± 0.1 1.5 ± 0.1 



1 sr 1 ) evaluated at the mid- bin 



"MeV 

b E? f J{E 1 ) (10~ 4 MeV cm- 
energy. 

c Total/7r°-decay/bremsstrahlung/inverse Compton. 
d The GALPROP galdef ID for this model is 54_5gXva rh7S which 
is available at the website http://galprop.stanford.edu 
e Statistical errors only. 
^Unidentified background. 



for 7-ray conversions in the front section of the LAT, and 
the sources are less well-localised spatially. In addition, 
the sky maps are well populated even at high latitudes 
and display spatial structure. The PSF broadening of the 
sources provides spatial structure and because the DGE 
model is fixed, more intensity is assigned to the source 
component to compensate in the fit. These effects lead 
to the systematic error in the source component being 
relatively larger than the isotropic at low energies and 
vice versa at high energies. Note, this applies for the 
high-latitude region from where the UIB component is 
derived, and also for the mid-latitude range for which 
we show the combined contribution by sources in Fig. [5] 
Because the uncertainties in the source and UIB compo- 
nents are not independent we have conservatively added 
their systematic uncertainties for the total intensity band 
shown in Fig. [2] 

The UIB component comprises the true extragalactic 
diffuse 7-ray emission, emission from unresolved Galac- 
tic and extragalactic sources, and residual particle back- 
grounds (CRs that pass the 7-ray classification analysis 
and 7-rays produced by CR interactions in the passive 
material outside the ACD) in the LAT data. In addi- 
tion, other relevant foreground components that are not 
completely modelled, such as emission from the solar disk 
and extended emission [22] ] and other potentially relevant 



"diffuse" sources [23j are included. Hence, the UIB com- 
ponent does not constitute a measurement of the extra- 
galactic diffuse emission. Furthermore, comparison with 
the EGRET estimate of the extragalactic diffuse emis- 
sion 24| is problematic due to the different DGE models 



used and analysis details that are beyond the scope of 
the current paper and will be addressed in a subsequent 
publication [251 ]. 

Discussion: The intensity scales of the LAT and 
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EGRET have been found to be different with the re- 
sult that the LAT-measured spectra are softer. In our 
early study of the Vela spectrum (26[, which was made 
using pre-launch IRFs, the difference was apparent al- 
ready above 1 GeV. Following on-orbit studies new IRFs 
have been developed to account for inefficiencies in the 
detection of 7-rays in the LAT due to pile-up and ac- 
cidental coincidence effects in the detector subsystems. 
The inefficiency increases at lower energies, with the re- 
sult that the IRFs used in the present analysis indicate 
greater intensities in the range below 1 GeV, with the 
magnitude of the effect ranging up to ~30% at 100 MeV. 
A forthcoming study of the Vela pulsar using the LAT 
one-year data with post-launch IRFs also shows a similar 
effect in the low-energy pulsed spectrum. So, the relative 
brightness of the diffuse emission measured by the LAT 
at low energies is unlikely to be due to increased residual 
background. Our confidence that the IRFs used in the 
present analysis accurately represent our knowledge of 
the instrument comes from detailed instrument simula- 
tions that were validated with beam tests of calibration 
units, and to post-launch refinements based on actual 
particle backgrounds. The systematic uncertainty on the 
effective area gives an energy dependent measure of our 
confidence in the IRFs used in the present analysis. 

As a consequence, the LAT-measured DGE spectrum 
averaged over all Galactic longitudes for the latitude 
range 10° < \b\ < 20° is systematically softer than the 
EGRET-measured spectrum. The spectral shape is com- 
patible with that of an a priori DGE model that is con- 
sistent with directly measured CR spectra. The excess 
emission above 1 GeV measured by EGRET is not seen 
by the LAT in this region of the sky. 

While the LAT spectral shape is consistent with the 
DGE model used in this paper, the overall model emis- 
sion is too low thus giving rise to a ^ 10 — 15% excess 
over the energy range 100 MeV to 10 GeV. However, the 
DGE model is based on pre- Fermi data and knowledge 
of the DGE. The difference between the model and data 
is of the same order as the uncertainty in the measured 
CR nuclei spectra at the relevant energies 27] . In ad- 
dition, other model parameters that can affect the 7-ray 
production rate (e.g., the conversion between CO line in- 
tensity and molecular hydrogen column density in the 
interstellar medium, Xqo) have not been modified in the 
present paper. Overall, the agreement between the LAT- 
measured spectrum and the model shows that the fun- 
damental processes are consistent with our data, thus 
providing a solid basis for future work understanding the 
DGE. 
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